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The losses, phase shifts and field distribution functions for the two lowest- 
order modes of interferometer-type maser resonators consisting of spherically 
curved mirrors with circular apertures are computed by solving a pair of in- 
tegral equations numerically on a digital computer. Solutions are obtained 
for the symmetric geometry of identically curved mirrors and for the half- 
symmetric geometry consisting of one plane and one curved mirror, with the 
radius of curvature of the mirrors as a variable parameter. The confocal or 
near-confocal configuration is shown to have good mode-selective properties 
in that the ratio of the loss of the second lowest-order (TEM W ) mode to that 
of the lowest-order (TEMoo) mode is the largest of the configurations con- 
sidered. The numerical results should be of interest to those concerned with 
the problem of mode selection in optical masers and with the design of single- 
mode maser s ivith relatively low gain. 

I. INTRODUCTION 

Interferometer-type resonators used for optical masers usually have a 
number of modal resonances falling under the gain profile of the active 
medium. Therefore optical maser oscillators generally can and often do 
oscillate in many modes, each mode having its own characteristic fre- 
quency and field pattern. Such a multimode, multifrequency output is 
undesirable for applications in communications and metrology. Many 
mode-selection schemes have been devised and tried, but most of them 
involve added complications and are beset by problems of stability. A 
simple solution was obtained by Gordon and White, 1 who built a stable 
single-frequency gas optical maser using a short, thin discharge tube. 
The resonator length was made short to reduce the number of longi- 
tudinal resonances and the mirror curvature was chosen so that only the 
lowest-order transverse mode had enough gain for oscillation. (The dif- 
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fraction losses of the higher-order modes were all greater than the gain 
of the tube.) Choosing the appropriate mirror curvatures for mode sup- 
pression is preferable to the commonly used method of aperturing the 
mirrors, because the former method does not restrict the amount of ac- 
tive material that can participate in maser action and therefore is capa- 
ble of producing greater output power. 

In order to select a pair of mirrors with the appropriate radii of curva- 
ture for single-mode operation, it is important to know accurately the 
diffraction losses of the modes as functions of the mirror size, spacing and 
curvature. Much work has been done on resonators with rectangular 
mirrors: Boyd and Gordon 2 and Boyd and Kogelnik 3 have derived ap- 
proximate formulas for estimating the losses; Fox and Li 4 have made 
numerical calculations on a digital computer by solving the appropriate 
integral equations by the method of successive approximations; Streifer 6 
recently made similar calculations using Schmidt expansion theory; 
Gloge 6 used perturbation techniques for his calculations. Unfortunately, 
these results are not very practical because the discharge tubes usually 
have circular cross sections and the diffraction losses for square mirrors 
and for circular mirrors are not simply related. Recently, Heurtley 7 has 
made some calculations on the modes and the eigenvalues of resonators 
with circular mirrors using Schmidt expansion theory. But published 
data on the diffraction losses suitable for practical design purposes are 
still lacking. It was with this purpose in mind that we computed the 
diffraction losses of the two lowest-order modes of the maser resonator 
with circular mirrors for various mirror curvatures. The results, which 
complement our earlier work, 48 - 9 should be of interest to those concerned 
with the problem of mode selection in optical masers. 

II. MATHEMATICAL FORMULATION 

The geometry of a maser resonator with circular mirrors is shown in 
Fig. 1. It is convenient to use the cylindrical coordinate system, with its 
axis coinciding with the axis of the resonator. The mirrors are spherically 
curved and have radii of curvature equal to Ri and R% . The mirror aper- 
tures are circular, their radii being a\ and a 2 . The separation between the 
mirrors along the axis is d. 

The integral equations for the modes of the resonator with mirrors of 
arbitrary curvature and shape are given in Ref. 4. For the present case, 
they are in the form 

t (1 V (1) (n , vO = r ( * K ™ fo >& ; r * 1 w) * (2) ^ ' ^ r *W r * (1 > 
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MIRROR I 
(RADIUS OF CURVATURE = R.) 



MIRROR 2 
(RADIUS OF CURVATURE = R 2 ) 



Fig. 1 — Geometry of the maser resonator consisting of spherically curved 
mirrors with circular mirror apertures. 

7 (2 V (2) fo , w) = P [ K a) fa ,<P2-,n, *,,) * (1) (ft , ^i) ndpidn (2) 

•'o •'o 
where K {1) and K i2) are equal and are given by 
K m (ri , vi ; r 2 , <p 2 ) = K«(n , <p 2 ; r, , *>,) 

= -^ exp < - 3 — [gn'i + 2 r 2 2 - 2rir 2 cos (pi - y 2 )] > • 

In the above equations, 7 (1) and y (2> are the eigenvalues associated with 
the eigenfunctions \l/ a) (ri,tp\) and ^ (2> (r 2 , v» 2 ), which are distribution 
functions of the reflected field at each mirror surface; g>i is equal to 
1 — (d/Ri) and g 2 is equal to 1 — (d/R 2 ) ; A; is 27r/X and X is the wave- 
length in the medium between the mirrors. Making use of the relation 10 

exp[jn(w/2 - P))Jn(xy) 

= ^— / exp j jU'y cos (a — 0) — na] } da 
•S7T •'o 

and integrating (1) and (2) with respect to <pi and <p 2 , respectively, it is 
seen that 



and 



* a >(ri,*i) = ^ (l) (n)e-'"*i, (n = integer) (5) 

^ (2) (r 2 ,ip->) = R» m fy*)*~ i "h (n = integer) (6) 



920 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1965 

satisfy (1) and (2), respectively. Thus the azimuthal variation is sinus- 
oidal in form. The radial functions fl„ (1) (rO and R n (2) (r 2 ) satisfy the re- 
duced integral equations 

y, t a) Rn W to Vn = r K n (r, , r 2 ) R n m (r 2 ) y/7, dr 2 (7) 

•'o 

y n (2) R™ ft) y/~r, = P K n (n , r 2 ) R» (1) (n) Vn dn (8) 



where 



(9) 



Kn ( ? , , r2 ) = ££. J n (k ™) Vw exp {- | W + «, 2 r 2 2 ) } 

and J„ is a Bessel function of the first kind and nth order. It is unneces- 
sary to show that the modes are orthogonal over their respective mirror 
surfaces, since this has been demonstrated for the general case of arbi- 
trary mirrors. 4 

III. COMPUTED RESULTS AND DISCUSSION 

An IBM 7094 computer was programmed to solve (7) and (8) for the 
two lowest-order modes (TEMoo and TEMio) and their eigenvalues using 
the method of successive approximations. 4 - 8 Solutions were obtained for 
the symmetric and half -symmetric geometries; the symmetric geometry 
consists of identical mirrors (ai = a 2 = a, gi = g% = g) and the half- 
symmetric geometry consists of one plane and one curved mirror 
(ai = a 2 = a, gi = 1.0). In all cases, fifty or more intervals were used 
for the numerical integration of (7) and (8). When the losses were low 
and the convergence was slow, as many as one hundred intervals were 
used. 

Figs. 2 and 3 show the relative field distributions of the lowest-order 
(TEMoo) and the next lowest-order (TEMio) modes for the symmetric 
case. Except for geometries close to plane-parallel or concentric (g w ±1), 
the relative amplitude distributions can be closely approximated by 
Gaussian-Laguerre functions, with the spot sizes of the TEMoo mode 
equal to (Kd/ir)^! — g 2 )~ l , as given by Boyd and Gordon. 2 Also, for 
geometries other than plane-parallel or concentric, the equiphase sur- 
faces almost coincide with the mirror surfaces. The agreement between 
the computed field distributions and those of the generalized confocal 
theory 2 ' 3 becomes closer with larger Fresnel numbers (N > 1). 

The power losses per transit of the two lowest-order modes are given 
in Figs. 4 and 5 for various values of g lying in the low-loss region 
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Fig. 2 — Relative field (list ril nit ions of the fundamental (TEMon) mode for 
the symmetric geometry (N - ar/\d = 1.0). 
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Fig. 3 — Relative field distributions of the TEMio mode for the symmetric 
geometry (N = a 2 /\d = 1.0). 
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Fig. 4 — Power loss per transit of the fundamental (TEMoo) mode for the 
symmetric geometry. 



(0 ^ | g | ^ I). 3,4 The curves for | g \ = 1 (plane-parallel or concentric) 
and g = (confocal) are the same as those given previously by Fox and 
Li. 8 Curves for other values of g lie between these two. The computed 
loss values for the confocal configuration are in perfect agreement with 
those of Slepian, 11 who has obtained various analytical expressions to ap- 
proximate the losses. For our interest, the most useful form of the ex- 
pression for the loss of the TEM nm mode obtained by Slepian is 



loss = 



27r(87rA0 



2»i+n+l -irN 
6/ 



r (m + l)T(m +n+l) 

constant 



_ (10) 



1 + 



2irN 



higher-order terms in 



2ttJV 
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Fig. 5 — Power loss per transit of the TEM 10 mode for the symmetry geometry. 



where n and m are the azimuthal and the radial mode numbers respec- 
tively, and T(.r) is the gamma function. The computed curves for the 
plane-parallel or the concentric configuration are in excellent agreement 
with those of Vainshtein, 12 who has obtained an approximate formula for 
the loss of the TEM„„, mode. It is of the form 



loss = 8v„, 



(11) 



[(M + 0) 2 + /3 2 ] 2 

where v nm in the rath zero of the Bessel function J n (x), /3 = 0.824 and 
M = (8tN)K Both formulas (10) and (11) fail for small Fresnel num- 
bers where losses arc appreciable. 

It is seen from Figs. 4 and 5 that the diffraction losses are very sensi- 
tive to changes in mirror curvature ; also, provided | g | is not very close 
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to unity, the losses decrease very rapidly as N increases. By choosing g 
and A r suitably, it is possible to have a strongly oscillating fundamental 
mode in an optical maser with all higher-order modes suppressed. For 
example, if a maser tube has a gain of five per cent and the Fresnel num- 
ber of the resonator is two, higher-order modes will not oscillate with 
g = 0.90. (Loss of TEM o mode ^1 per cent and loss of TEMio mode 
£^5.2 per cent). 

While the loss of each mode is given by the magnitude of its eigenvalue, 
its resonant frequency is determined by the phase of its eigenvalue, 
which is equal to the phase shift (relative to the geometrical phase shift) 
per transit for the mode. The phase shifts for the two lowest-order modes 
are plotted in Figs. 6 and 7. The curves shown are for positive g only; 



100 
80 

60 
40 



1.0 
0.8 



0.4 

























































5ii ^; 












-g-o 


















^^s 


^ 










0.5 






























0.9 










































































\i.o 






















































































































TEM 00 MODE 







































































































































































2 4 6 8 10 

N = a 2 /Xd 



Fig. 6 — Phase shift per transit (leading relative to the geometrical phase shift) 
oi the fundamental (TEM W ) mode for the symmetric geometry. 
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Fig. 7 — Phase shift per transit (leading relative to the geometrical phase 
shift) of the TEM™ mode for the symmetric geometry. 

the phase shift for negative g is equal 4 to 180 degrees minus that for posi- 
tive g. The horizontal portions of the curves can be calculated from the 
theory of Boyd and Kogelnik 3 as 



phase shift = (2m + n + 1) arc cos VWa 

= (2m + n + 1) arc cos g, for g\ = gi • 



(12) 



The problem of mode discrimination in an optical maser can be ap- 
proached in two different ways which depend on the available gain of 
the active medium. If the gain is large, the relative magnitudes of the 
eigenvalues of the different modes are of importance. 131416 On the other 
hand, if the gain is small, such as in 6328A He-Ne masers, the relative 
losses are important. 9 The ratio of the loss of the TEMio mode to the loss 
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Fig. 9 — Average power loss per transit of the fundamental (TEMoo) mode for 
the half-symmetric geometry. 
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Fig. 10 — Average power loss per transit of the TEMio mode for the half- 
symmetric geometry. 

of the TEMoo mode, which is a measure of mode selectivity for the low- 
gain case, is plotted in Fig. 8 as functions of the Fresnel number. The 
dotted curves are contours of constant loss for the TEMoo mode. It is 
clear from this plot that the confocal or near-conf ocal geometry possesses 
good mode-selective properties, and that the plane-parallel or the con- 
centric geometry is rather poor. However, if an aperture is placed in the 
midplane of a concentric resonator, its mode-selective properties have 
been shown to improve and to approach the confocal resonator. 9 

One of the commonly used resonator configurations is the half-sym- 
metric geometry consisting of one plane and one curved mirror. If the 
spot size of the mode pattern at the plane mirror is very much smaller 
than the effective mirror aperture, such as in the case of the half-con- 
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centric or near half-concentric configuration, the half -symmetric resona- 
tor can be regarded as equivalent to a symmetric resonator of twice the 
length. However, if the mirror configuration is such that it lies between 
the half-confocal and the plane-parallel (g t =1.0 and 0.5 ^ g? ^ 1.0) 
and the effective mirror apertures are the same, the equivalence of half- 
symmetric and symmetric resonators is not valid and the modes and their 
eigenvalues must be recomputed. Computations were carried out for 
0i=l.O and 02 = 0.5, 0.75, 0.9 and 0.95 corresponding to = 0, 0.5, 
0.8 and 0.9 of the symmetric case. Since the mode pattern at the plane 
mirror is now different from that at the curved mirror, the losses and the 
phase shifts incurred at the mirrors are also different. As before, 4 we de- 
fine the average loss per transit as 1 — | 7 (1 V 2) | and the average phase 
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Fig. 11 — Average phase shift per transit (leading relative to the geometrical 
phase shift) of the fundamental (TEMw) mode for the half-symmetric geometry. 
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shift per transit as (phase of 7 (1) + phase of y (2) )/2. The computed 
average losses and average phase shifts per transit for the two lowest- 
order modes are given in Figs. 9 to 12, which are to be compared with 
those of the symmetric geometry shown in Figs. 3 to 6. The ratio of the 
loss of the TEMio mode to the loss of the TEMoo mode for the half -sym- 
metric resonator is plotted in Fig. 13. Comparing Figs. 8 and 13, it is 
seen that the mode-selective properties of the half-symmetric resonator 
are very similar to those of the corresponding symmetric resonator. 

Although the computed results are for equal mirror apertures (ai = a 2 ) , 
they are also applicable to certain "equivalent" geometries with un- 
equal mirror apertures. The equivalence relations are discussed and 
given in a recent paper by Gordon and Kogelnik. 16 
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Fig. 12 — Average phase shift per transit (leading relative to the geometrical 
phase shift) of the TEMio mode for the half -symmetric geometry. 
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Fig. 13 — Ratio of the average losses per transit of the two lowest-order modes 
for the half-symmetric geometry. The dotted curves are contours of constant 
average loss for the TEMoo mode. 

IV. CONCLUSIONS 

The diffraction losses of the modes of interferometer-type maser 
resonators can be utilized for mode discrimination. Results of the com- 
putation show that the confocal or near-confocal resonator has good 
mode-selective properties; but its mode volume is small compared with 
those of the plane-parallel or concentric resonator of the same Fresnel 
number. 2 • l7,18 Since it is desirable to utilize as much of the active mate- 
rial as possible in order to obtain the maximum output power from a 
maser, a geometry that is as close to plane-parallel as is consistent with 
the requirements of mode discrimination and mechanical stability is 
to be preferred. In practice, the perturbing effects of the maser tube 
wall, the mirror irregularities, the nonlinearity of the active medium, 
etc. will modify the idealized modes and their losses and therefore should 
be taken into consideration. 
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